Magnetic Field Effect in Josephson Tunneling between d-Wave Superconductors 

Xin-Zhong Yan^'^ and Chia-Ren Hu^ 
^Department of Physics, Texas A&M University, College Station, TX 77843-4242 
^Institute of Physics, Chinese Academy of Sciences, P. O. Box 603, Beijing 100080, China 
Q^". (February 1, 2008) 

0^ 



The magnetic field effect in the Josephson tunneling between two d-wave superconductors are 
investigated. When the crystal orientation of one (or each) superconductor relative to the interface 
normal is such that midgap states exist at the interface, there is a component of the tunneling 
1-^ ' current due to the midgap states. For a junction with a flat {001}|{110} or {100}|{110} interface, 

this component is the predominant contribution to the current. The predicted current-field de- 
pendence differs entirely from the conventional Fraunhofer pattern, in agreement with a published 
^O ' measurement. This is because, apart from the Fraunhofer factor, the critical current depends on 

Cn , the magnetic field B through the current density also which is a linear function of B for weak B. 
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Whether the magnetic-field dependence of the Joseph- 
son tunnehng current between two high-Tc superconduc- 
tors (HTSC's) should show the well-known Fraunhofer 
pattern seen in junctions made of conventional supercon- 
ductors is still an unsettled question. There seem to be no 
theoretical discussions sheding light on this subject so far. 
On the other hand, many measurements have been de- 
voted to the investigation of this problem [1]. However, it 
is not easy to fabricate high-quality junctions of HTSC's 
with good flat interfaces. Since the pairing order pa- 
rameter in HTSC's very likely has predominately d-wave 
symmetry, the behavior of the Josephson tunneling cur- 
rent is expected to depend strongly on the interface ori- 
entations. Most of the Josephson junctions exhibit ran- 
dom patterns in their current-field relations most-likely 
because their interfaces are made of many facets with 
a variety of orientations [1,2]. Recently, Ishimaru et al. 
have successfully fabricated junctions with atomically fiat 
{001}|{110} interfaces [3]. Their measurements showed a 
regular dependence of the Josephson critical current on 
the applied magnetic field, which is fundamentally dif- 
ferent from the Fraunhofer pattern. Such a seemingly 
anomalous behavior of the Josephson current has not yet 
been explained in an acceptable theory. 

In this letter, we present a calculation of the Joseph- 
son tunneling current between two d-wave superconduc- 
tors (DWSC's) with a {001}|{110} interface, as is inves- 
tigated in that experiment. Our calculation shows that 
the local Josephson-current density depends linearly on 
the magnetic field when the field is weak. Therefore, the 
current-field relation in such a Josephson tunneling junc- 
tion is predicted to be different from the usual Fraunhofer 
pattern. This behavior stems from the existence of the 
midgap states (MS's) at the junction interface [4-8]. 

It has been noted that the MS's play a very impor- 
tant role in Josephson tunneling processes between two 
DWSC's [9] (and also in quasi-particle tunneling injunc- 
tions involving DWSC's which does not concern us here). 
Consider a SC with a {110} surface at x := 0, with the 
X axis being defined perpendicular to this surface. Let a 
magnetic field B be applied along the c-axis, so that there 
is a vector potential A along [110], which is designated 
as the y axis. In the gauge in which the order parameter 
is real, the vector potential vanishes in the interior of the 
SC. We have A{r) ~ -BAexp(-a;/A), where A ~ 1500l 
is the magnetic-field penetration depth in the ab plane. 
The perturbation to the motions of particles is described 
by the Hamiltonian, 



H' ^-- I d^rj{r)-A{r). 



(1) 



where c is the speed of light, and j(r) the current den- 
sity. The particles in the continuum states should only be 
perturbed very weakly by this magnetic field since their 
wave functions extend far away from the surface where 



S ~ 0. In contrast, the perturbation to the MS's is sig- 
nificant because the MS's are localized within roughly a 
coherence length (^ ^ 15A) from the surface, which is 
much smaller than A. From the Bogliubov-de Genncs 
equations for quasi-particles which are coherent super- 
positions of electrons and holes, the eigen-energies of the 
perturbed MS's can be obtained as 
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where k is the wave number along y, and we have ne- 
glected the Zeeman energy. [10] This energy expression 
is an odd function of k. Thus at sufficiently low temper- 
atures, an asymmetric distribution of quasi-particles in 
the MS's can be realized. 

In the Josephson junction under consideration, the 
wave functions of the MS's are distributed predominately 
in the side with a {110} surface, with only a small ampli- 
tude tunnelling into the other side (assuming high tun- 
neling barrier at the interface). These states are Andreev 
bound states localized within roughly a coherence length 
from either side of the interface. It is well known that 
quasi-particles in such bound states can transport electric 
current between two superconductors. In the elementary 
charge-transport process in this problem, the direction 
of the tunneling current mediated by such a state will 
depend on the sign of its momentum k along y, since this 
tunneling current is proportional to the wave-function 
product u^Vk which is proportional to the order param- 
eter of this SC at its surface, the sign of which depends 
on the sign of fc. At zero magnetic field two states of 
momentum k and — fc, respectively, have the same en- 
ergy E{k) — E{—k). They arc thus equally occupied, so 
that their contributions to the Josephson current cancel 
each other. Under a small non-vanishing magnetic field, 
the energies of two MS's of such a pair of fc values as 
given by Eq. (2) arc opposite in sign. Their occupa- 
tions at low temperatures will no longer be equal. With 
cancellation no longer complete, a net current can now 
be transferred across the Junction via this pair of MS's 
which is proportional to the difference of their occupation 
./(£fc) - ./(e-fe) = tanh(e_fc/2fcBT), where /(ej.) is the 
Fermi distribution function at temperature T and fcs is 
the Boltzniann constant. If the magnetic field applied in 
the tunneling experiment is so weak that \ek\ <C 2kBT — 
roughly, the ratio |efe|/2fcBT is less than 0.05-B/T where 
B and T are in units of Gauss and Kelvin, respectively, 
— the current density J{B) is then a linear function of 
B. The total critical current Ic{B) of a Josephson junc- 
tion is determined by this current density times a factor 
describing a Fraunhofer-like pattern. As the final result, 
at weak magnetic field, Ic{B) is totally different from the 
usual Fraunhofer pattern. 

To quantitatively study this problem, we use the 
tunneling-Hamiltonian approach in a tight-binding 
model. For the sake of description, the SC side with 



a {110} ({001}) surface is referred to as the ridgt (left) 
side of the junction. The bridge width of the junction is 
denoted as L. As is the situaton in the experiment, we 
assume that the barrier thickness and the length of the 
left SC in the direction perpendicular to the interface are 
both much less than A. Then the effects of the magnetic 
field in this side of the junction and in the barrier region 
can be neglected. Note that the external magnetic field 
Bext applied to the junction can induce a paramagnetic 
field Bpara to the right of the surface region occupied by 
the MS's in the right SC, so that the total magnetic field 
applied to the right SC is practically B — Bext + Bpara, 
which is then screened in a distance ^ A in the right SC. 
According to the tight-binding model adopted here, the 
electrons are supposed to tunnel from the surface atoms 
of one SC to the neighboring atoms on the surface of the 
other, through the interface. The Hamiltonian describing 
this tunneling process is given by 



^T-5](r^,dLc,o + ^-c.) 



(3) 



where dj^, (cj^) is the creation (annihilation) operator 
for electrons with spin-a in the right (left) SC, the i 
and j summations run over the surface sites on the two 
sides of the junction, respectively. (In the actual cal- 
culation, each sum is over one infinite chain of atomic 
sites only, since for the right side we need only consider 
one Cu02 plane with one linear {11} edge, which practi- 
cally only interacts with one linear {10} chain of atoms 
on the left side.) The tunneling matrix element is then 
Tij = jTijI exp[i(0o/2 -l-fcoyi)], with k^ = eBX/Uc and (f)a 
is the initial phase difference between the two SC's. The 
magnitude of the tunneling matrix element \Tij\ is con- 
sidered as a function of the distance \yi — yj\ (where yi 
and yj are respectively the coordinates of the sites i and 
j along the interface). Here, we suppose a simple func- 
tional form |Ty| = T(j/j - y^) = Toexp(-(yi - yjf /rl) 
with the parameter tq describing the effective tunneling 
range. (The precise functional form is not important as 
long as it is short-ranged. The chosen form is for con- 
venience only.) By the standard perturbation treatment, 
the Josephson tunneling current can be expressed as 



/ = 2elr[v<^r{uj)\^^Q 



(4) 



with <I>r('-^) the analytical continuation iuon (= i2mTkBT) 
-^ Lu + iO^ of the correlation function $,.(ia;„). In terms 
of the one-particle Green functions. 
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(5) 
where fte = (2£ + l)TrkBT is the Matsubara frequency, 
and Fln{Vtt) and F- , (fi^) are the Fourier transforms 
of the imaginary-time (t) anomalous Green functions 
Fkij) = - < T,dl^{r)d\^{Q) > and F^^,{t) = - < 



Tt-cj^{t)cj' i{0) >, respectively. Along the interface of 
the junction, the anomalous Green functions depend on 
the relative distance of two y-coordinates only. Ac- 
cordingly, we factorize TijTi'ji to \TijTi>j>\exp[iko(yi> — 
yi)]exp[i((/)o -|- 2fcoJ/i)]- The factor exp[ifco(j/i' - yt)] may 
be approximated by unity, since it leads to a negligi- 
ble correction to a momentum conservation relation ob- 
tained later when Fourier transform is employed, (fco is 
estimated to be ~228 cm^^ at _B = IG, whereas the 
important momenta are all of order kp.) The summa- 
tion over i which can be performed immediately then 
produces the factor describing a Fraunhofer-like pattern: 
jo{koL) = s[n{koL)/koL. 

To illustrate the physics, we neglect the space varia- 
tion of the order parameter near the interface since the 
existence of MS's is invariant with respect to such vari- 
ations. With this simplification, the wave functions of 
the MS's and also the Green's functions can be obtained 
analytically within the tight-binding model defined in a 
square lattice [8,11]. In momentum space, the anomalous 
Green's function F' is given by 



F\k,ni) = -a 
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where a is the lattice constant, Ag — 2A[cos{qxa) — 
cos(fca)] and Eq [with q = {q^, k)] are the order parameter 
and the quasi-particle energy, respectively [8,10]. For the 
problem under consideration, we need only consider the 
current transferred via the MS's. Therefore, the anoma- 
lous Green function F'' is simply: 
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where u^ and v^. are the surface values of the wave 
functions of the MS with momentum k. In the lat- 
tice model, the eigen-energy of a MS is given by Sk = 
— {2^/2ateB\/hc) sm{ka/V2), and t is the hopping en- 
ergy of electrons between nearest-neighbor sites. It is 
worth noting that Vk = WfeSgnfc for these MS's, which 
means that the tunneling-current direction depends on 
the sign of k as mentioned above. Also, due to the phase 
factor i, the phase dependence of the Josephson current 
is cos (j)o rather than the usual sin 0o • Substituting these 
results into Eq. (4), we obtain 



/ = J{B)jo{koL)cos(l)o, 



(8) 
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where Ns is the number of total lattice sites on the {110} 
surface of the right SC, T{k) is the Fourier transforma- 
tion of T{y), q = {qx, k + Qn) with Qn ~ 2n7r/-\/2a, and 
g{ek,Eq) = [Unhiek/2kBT)-te.nhiEj2kBT)]/iek-Eg). 
Apart from the factor jo(fcoi), the Josephson current de- 
pends on B through J{B) also. At weak B, J{B) ex B. 
This is an unusual feature of the Josephson current be- 
tween two DWSC's when MS's are present at the inter- 
face and give the dominant contribution to the current. 
In Fig. 1, we show the numerical results for J{B) nor- 
malized by Jo = eNsT^A'^/t^ at T = 4.2K for a junction 
with a {001}|{110} interface. The long-dashed, solid, 
and short-dashed curves correspond to the parameters 
qq = Tq^ = 1.0a~^, 0.94a~^, and 0.9a~^, respectively. 
Obviously, J{B) ex B at weak B. The slope of J{B) 
&t B = depends sensitively on the parameter tq. For 
a narrow region in the vicinity of qo = 0.94a~^, J{B) 
changes sign at a certain value of B of the order of Bq. 
This is because right at the interface, the pairing am- 
plitude in the left SC is proportional to Acos(fco). Its 
sign and therefore the sign of the corresponding con- 
tribution to the tunneling current density depend on 
the magnitude of the momentum k of the quasi-particle: 

cos(A;a) < for ka > 7r/2. Owing to the _B-dependence 
of the energy of a MS, \ek\ oc Bsm{ka/V2), changing 
magnetic field is equivalent to modulating the distribu- 
tion of the quasi-particles in the k space. Similarly, the 
relative weights of positive and negative currents are also 
controled by rg. These effects result in the delicate de- 
pendence of J{B) on B and rg. 

There is some subtlety in the factor jo{koL) as well: In 
an ideal system, there exists a spontaneous current near 
the {110} surface of the right SC [10,12,13], so a B^xt- 
independent spontaneous paramagnetic field Bpara = 
Bo = hc/2e}? (which is ^146 Gauss as estimated by 
using A ~ 1500 A) is produced. The total magnetic 
field B ~ Bpxt + Bq in the right SC is nonzero even at 
Bext = 0. Therefore, the factor jo(fcoi) alone can already 
cause the Ic-B^xt curve to not have the usual Fraunhofcr 
pattern. In particular, there may be no main peak at 
Bext = as in the usual Fraunhofcr pattern. Generally, 
there is no unique pattern if this spontaneous field exists: 
At Bext = 0, there may be a dip or a cusp in Ic{Bext) 
depending on the geometric ratio L/X of the junction. 
However, the existence of this spontaneous magnetic field 
depends sensitively on the quality of the juction includ- 
ing the flatness of the interface and the impurity content 
of the SC which has the MS's. We do not believe that 
most of the HTSC's can have this spontaneous current. 
Even for a perfectly flat {110} interface, the existence of 
impurity scatterings is unavoidable. In such a non-ideal 
system, the particles in the MS's have finite lifetimes 
due to their scatterings off the impurities. Therefore, 
the spontaneous surface current and thereby the sponta- 
neous magnetic field cannot exist. The induced param- 



agnetic field is no longer a constant but depends linearly 
on the external field, Bpara oc Bext, for small field, when 
the magnitude of the maximum shift of the MS's is still 
much smaller than the width of the midgap-state peak 
in the density of states. Consequently, the dependence 
of Jo(fco-^) on the external field in the weak field regime 
is essentially the conventional Fraunhofcr pattern except 
for a scale factor (< 1) on B^xt- 

Shown in Fig. 2 is the result for Ic{B) normalized by 
Jo for the junction with L/X — 35 and qg = 0.94a~^. To 
compare with the experimental measurements, we need 
to take into account a larger enhancement factor in the 
effective magnetic field, which is mainly due to the flux 
focusing effect in thin-film grain boundary junctions [14]. 
Due to this effect, the effective field B can be an order of 
magnitude larger than the external field B^xt- Under this 
consideration, the experimental result of Ishimaru et al. 
[3] can be explained by the present theory, as is shown 
in Fig. 2. According to the present theory, the deep 
dip in the measured I{B) curve at zero magnetic field is 
a strong indication that the pairing has predominately 
d-wave symmetry in the Cu02 planes of YBCO. 

If the SC's have the d+s symmetry, then the tunneling 
current can also be transferred through the continuous 
states. In this case, the Josephson critical current Ic{B) 
does not vanish at _B = 0. There may appear a dip or 
a peak in the I — B curve at _B = depending on the 
magnitude of the s-wave component. The fact that the 
observed /c(0) is small with no peak there should imply 
that the s-wave component is small, though possibly not 
zero in YBCO. The predicted features of the Ic — B curve 
can also be observed in junctions with {100}|{110} inter- 
faces if only they can be fabricated with flat interfaces. 
The physics is essentially the same as for the case of a 
{001}|{110} interface. 

In summary, we have investigated the magnetic-field 
dependence of the Josephson current in a junction made 
of two DWSC's with a {001}|{110} interface. In such a 
junction the Josephson current is mainly mediated by the 
MS's on the side with a {110} surface. Under a magnetic 
field, the eigen-energies of these states shift in direct re- 
lation to their momenta along the interface, leading to 
a skewed occupation of these states, which implies a B- 
dependent current-density factor J{B) which is linear in 
B for small B. When the particles in the MS's have a fi- 
nite life-time, the effective field B acting on the side with 
a {110} surface should be a linear function of the external 
field when the field is weak. Thus the usual Fraunhofcr 
factor is not modified qualitatively, but the critical cur- 
rent Ic depends on B through a current-density factor 
J{B) also which varies linearly with B at weak B. Thus 
the Ic — Bext curve differs entirely from the conventional 
Fraunhofcr pattern, especially near zero field. 

This work is supported by the Texas Higher Education 
Coordinating Board under the grant No. 1997-010366- 
029, and by the Texas Center for Superconductivity at 



the University of Houston. 



[1] N. G. Chew, et al, Appl. Phys. Lett. 60, 1516 1992); 

B. V. Vuchic, et al, J. Appl. Phys. 77, 2591 (1995); J. 

Mannhart, et al, Phys. Rev. Lett. 77, 2782 (1996); H. 

Hilgenkamp, et al, Phys. Rev. B 53, 14586 (1996). 
[2] R. G. Mints and V. G. Kogan, Phys. Rev. B 55, R8682 

(1997). 
[3] Y. Ishimaru, et al, Phys. Rev. B 55, 11851 (1997). 
[4] C.-R. Hu, Phys. Rev. Lett. 72, 1526 (1994). 
[5] J. Yang and C.-R. Hu, Phys. Rev. B 50, 16766 (1994). 
[6] L. J. Buchholtz, M. Palurnbo, D. Rainer, and J. A. Sauls, 

J. Low Temp. Phys. 101, 1079 (1995); 101, 1099 (1995). 
[7] M. Matsumoto and H. Shiba, J. Phys. Soc. Jpn. 64, 

1703 (1995); 64, 3384 (1995); 64, 4867 (1995); 65, 2194 

(1996). 
[8] X.-Z. Yan and H. lyetomi, Phys. Rev. B 57, 7944 (1998). 
[9] Y. Tanaka and S. Kashiwaya, Phys. Rev. B 53, R 11957 

(1996); and 56, 892 (1997); Yu. Barash, et al, Phys. 

Rev. Lett. 77, 4070 (1996); M. P. Samanta and S. Datta, 

Phys. Rev. B 55, R8689 (1997); R. A. Riedel and P. F. 

Bagwell, Phys. Rev. B 57 6084 (1998). 



[10] C.-R. Hu and X.-Z. Yan, cond- mat/9901334. This work 



has compared the sizes of the Zeeman and orbital shifts 

of the midgap-state energy. 
[11] X.-Z. Yan, Phys. Rev. B 56, 8374 (1997). 
[12] S. Higashitani, J. Phys. Soc. Jpn. 66, 2556 (1997). 
[13] M. Fogelstrom, D. Rainer, and J. A. Sauls, Phys. Rev. 

Lett. 79, 281 (1997). 
[14] P. A. Rosenthal et al, Appl. Phys. Lett 59, 3482 (1991). 

Figure Captions 

Fig. 1 The calculated current density J{B) of Joseph- 
son tunneling between two d-wave superconductors 
with a {001}|{110} interface as a function of the 
magnetic field B at T — 4:.2K. The long-dashed, 
solid, and short-dashed lines correspond to the pa- 
rameter Qo = r^^ — 1.0a~^, 0.94a~^, and 0.9a~^, 
respectively. The magnetic field is normalized by 
Bo = hc/2eX^. 

Fig. 2 The calculated critical Josephson current Ic{B) as 
a function of the magnetic field B for the junction 
of two d-wave superconductors with a {001}|{110} 
interface at T = 4.2^ and qo = OMa'^. The 
ratio between the bridge width and magnetic field 
penetration depth (in the ab plane) is L/X = 35. 
The black dots are data from Ref. 3. 
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